Experiments on Critical Phenomena in a Noisy Exit Problem 
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We consider noise-driven exit from a domain of attraction in a two-dimensional bistable system 
lacking detailed balance. Through analog and digital stochastic simulations, we find a theoretically 
predicted bifurcation of the most probable exit path as the parameters of the system are changed, 
and a corresponding nonanalyticity of the activation energy. We also investigate the extent to which 
the bifurcation is related to the local breaking of time-reversal invariance. 



PACS numbers; 05.40. -f-j, 02.50.-r, 05.20.-y 

Noise-induced motion away from a locally stable state, 
in a system far from thermal equilibrium, arises in di- 
verse scientific contexts, e.g., glasses [Q, arrays of Joseph- 
son junctions stochastically modeled computer net- 
works 1^, stochastic resonance [Q, and stochastic ratch- 
ets 1^. Because these systems in general lack detailed 
balance, progress in understanding this phenomenon has 
been slower than in thermal equilibrium systems. In par- 
ticular, there exist no simple or general relations from 
which the rate of noise-induced transitions between sta- 
ble states can be obtained. 

Recently, substantial progress on the nonequilibrium 
case has been achieved in the limit of weak noise, us- 
ing path integral or equivalent Hamiltonian formulations 
Fluctuational motion of the system can then be 
characterised by the pattern of optimal (i.e. most prob- 
able) fluctuational trajectories. An optimal trajectory 
is one along which a system moves, with overwhelming 
probability, when it fluctuates away from a stable state 
toward a specified remote state. These are rare events 
but, when they occur, they do so in an almost determin- 
istic way: e.g. escape from a domain of attraction typi- 
cally follows a unique trajectory . The properties of this 
most probable exit path (MPEP) determine the weak- 
noise behaviour of the mean first passage time (MFPT) . 

In recent years, it has been realized that in nonequilib- 
rium systems, the pattern of optimal fluctuational trajec- 
tories may contain "focusing singularities" [p|jic|]. Their 
effect on exit phenomena was considered by Maier and 
Stein |l^,|l^ who showed that, for a symmetric double 
well system (see below), the MPEP bifurcates when the 
model parameters are changed in such a way that a fo- 
cusing singularity appears along it. That is, the MPEP 
ceases to be unique. This bifurcation breaks the sym- 
metry of the model, and is accompanied by a nonan- 
alyticity of the activation energy for inter-well transi- 
tions: it is analogous to a second-order phase transition 
in a condensed matter system This analogy throws 
new light on, e.g., exit bifurcation phenomena in systems 
driven by coloured noise |Q| . 



Many of these theoretical ideas, although important, 
remain untested experimentally or numerically. In this 
Letter we use an analog experiment and numerical sim- 
ulations to demonstrate the predicted bifurcation of the 
MPEP, and the corresponding nonanalytic behavior of 
the activation energy and related quantities. We inves- 
tigate the nature of the broken symmetry in detail, and 
show how bifurcation is accompanied by a loss of time- 
reversal invariance along the MPEP. 

We investigate the motion of an overdamped particle 
in the two-dimensional drift field first proposed in [^: 
u(a;, y) — {x — — axy^ , —y — x'^y), where a is a param- 
eter. It has point attractors at (±1, 0) and a saddle point 
at (0,0). If the particle is subject to additive isotropic 
white noise f(i) = {fx,fy), its position {x,y) will satisfy 
the coupled Langevin equations 

x^x-x^ - axy"^ + /^(i), 

y = -y - x^y + fy{t), (1) 
(/^W) =0, {f,{s)f,{t))=D6,,5{s^t). 

Since u is not a gradient field (unless a = 1), the dy- 
namics will not satisfy detailed balance. The Fokker- 
Planck equation for the particle's probability density 
p ~ p{x, y, t) will be 

p={D/2)V^p-W-{pu). (2) 

In the weak- noise limit, escape of the particle from 
the domain of attraction of either fixed point (xs,0) — 
(±1,0) is governed by the the slowest-decaying nonsta- 
tionary eigenmode of the Fokker-Planck operator [p^ , 
whose eigenvalue Ai becomes exponentially small as 
D ^0. In this limit the MFPT (toxit) is well ap- 
proximated by Xi^. The slowest-decaying eigenmode is 
called the quasistationary probability density; we denote 
it by pi. It may be approximated in a WKB-like fash- 
ion |Jn|,|ll, i.e., 

pi{x,y)^Kix,y)e^pi-Wix,y)/D), D ^ 0. (3) 
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Here W{x, y) may be viewed as a classical action at zero 
energy, since it turns out to satisfy an eikonal (Hamilton- 
Jacobi) equation of the form i/(x, VW) — 0, where 
i7(x, p) = ^p^ + u(x) • p is a so-called Wcntzel-Freidlin 
Hamiltonian |Q. The optimal fluctuational trajecto- 
ries are projections onto coordinate space of the zero- 
energy classical trajectories determined by this Hamilto- 
nian. These lie on the three-dimensional energy surface 
specified by H — 0, embedded in the four-dimensional 
phase space with coordinates {x,y,px,Py)- In general, 
the computation of Vl^(a;, y) requires a minimisation over 
the set of zero-energy trajectories starting at [xg , 0) and 
terminating at {x,y). Moreover the MPEP, for the do- 
main of attraction of [xg, 0), is the zero-energy trajectory 
of least action which extends from {xs , 0) to the sad- 
dle (0,0). The MPEP action SW = VF(0,0) - W{xs,0) 
governs the weak-noise behavior of the MFPT. To lead- 
ing order it is of the activation type, i.e., 

(ioxit) - const X e^^/", D ^ 0. (4) 

So SW is interpreted as an activation energy. The prefac- 
tor 'const' is determined by the function K[x,y) of (||). 

When a = 1, the dynamics of the particle satisfy de- 
tailed balance, and the pattern of optimal trajectories 
emanating from (x^jO) contains no singularities. It was 
found earlier ]rT| , p^ that, as a is increased, the first fo- 
cusing singularity on the MPEP (initially lying along the 
X-axis) appears when a = ac = 4. It signals the ap- 
pearance of a transverse 'soft mode,' or instability, which 
causes the MPEP to bifurcate. Its physical origin is 
clear: as a is increased, the drift toward (a;s,0) 'soft- 
ens' away from the a;-axis, which eventually causes the 
on-axis MPEP to split. The two new MPEP's move off- 
axis, causing the activation energy (previously constant) 
to start decreasing. So the activation energy as a func- 
tion of a is nonanalytic at a — ac- If a is increased 
substantially beyond ac, further bifurcations of the on- 
axis zero-energy classical trajectory occur when a equals 
ckc"'^ = [j + 1)^, where j is the number of the bifurcation. 
But the oscillatory trajectories arising from such bifur- 
cations are believed to be unphysical, since the on-axis 
trajectory is no longer the MPEP. (Cf. p.) 

To test these theoretical predictions, and to seek fur- 
ther insight into the nature of the broken symmetry, we 
have built an analog electronic model of the system (|l|) 
using standard techniques [Tsl] . We drive it with zero- 
mean quasi-white Gaussian noise from a noise-generator, 
digitize the response x(t), y(t), and analyse it with a digi- 
tal data processor. Transition probabilities are measured 
by a standard level-crossing technique. Experimental in- 
vestigations of the optimal fluctuational trajectories are 
based on measurements of the prehistory probability dis- 
tribution 1 16 1^. This method was recently extended 
to include analysis of relaxational trajectories and thus 
to investigate directly the presence or absence of time- 
reversal symmetry and detailed balance pSl . 




FIG. 1. (a) The mean first passage time as a function of 
inverse noise intensity 1/ D for q = 1, from analog experiment 
(bars), numerical simulation (circles) and calculation (solid 
line), (b) The inter- well activation barrier 5W , as a function 
of a, from analog experiment (crosses), numerical simulation 
(squares) and theory [ p4]Jl^ (full line) . The dashed curve and 
circle data represent the MFPT prefactor from calculation 
and numerical simulation respectively. 



We have also carried out a complementary digital sim- 
ulation of (|^) using the algorithm of , with particular 
attention paid to the design of the noise-generator on ac- 
count of the long simulation times. Transition probabili- 
ties were measured using a well-to-well method, and the 
analysis of the data to extract the optimal fluctuational 
and relaxational trajectories was based on a method sim- 
ilar to that used in the analog experiments. 

Some activation energy results are shown in Fig. 1. 
Part (a) plots the MFPT (toxit) as a function of in- 
verse noise intensity 1/D for the special case a = 1. 
In this case the drift field is the gradient of the potential 
U{x, y) = (y2(l + a;2) - x'^ + x'^/2) /2, and W can be ob- 
tained exactly {W — 2U). The analog and digital results 
are in good agreement, and demonstrate that the noise 
dependence of the MFPT is indeed of the activation type 
predicted by the theory. Activation energies determined 
from the slopes of a series of plots like those in Fig. 1(a), 
yielded the results shown in Fig. 1(b), where they are 
compared with theoretical values of 6W determined from 
the true (least action) MPEP or MPEP's [|l]jl2|. At the 
predicted critical value ac = 4, marked changes in both 
the activation energy and MFPT prefactor (which 0) 
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FIG. 2. Measured positions of the ridges (first moments) 
for a = 6.67 of (a) the fluctuational part (filled circles) and 
(b) the relaxational part (pluses) of the escape trajectories. 
Nearby theoretical trajectories are shown by lines. The insert 
shows the exit location distribution along the y-axis. 

are evident: theory predicts that the activation energy bi- 
furcates here into two values, corresponding to paths on 
and off the x-axis, of which only the latter (lower action) 
path is expected to be physically meaningful. The depen- 
dence of the activation energy on a near the second crit- 

(2) 

ical value al ' = 9 is smooth, in agreement with the pre- 
diction that higher bifurcations correspond to folding of 
a nonphysical sheet of the 'action surface' W = W{x, y), 
and are not observable 

Interestingly, the transition shown in Fig. 1(b) resem- 
bles the bifurcation of the activation energy in an over- 
damped oscillator driven by coloured noise 0. This 
suggests that the WKB analysis |jll|,|l^ of (|l|) may pro- 
vide physical and topological insight into the correspond- 
ing transition phenomena in systems driven by quasi- 
monochromatic noise. 

To verify experimentally the expected relationship be- 
tween the bifurcation of the MPEP and the bifurcation of 
the activation energy, we have measured two-dimensional 
prehistory probability distributions of fluctuational 



trajectories bringing the system into the vicinity of the 
separatrix between the two wells (the y-axis) . In the limit 
of low noise intensity, the maxima of the corresponding 
distributions trace out optimal trajectories |l^Jl^. The 
positions of these maxima are compared to the calculated 
MPEP's for a = 6.67 in Fig. 2(a). It is clear that the 
typical fluctuational path corresponding to escape from 
the domain of attraction of {xs , 0) follows very closely 
one of the predicted MPEP's. 

To seek further experimental insight into the charac- 
ter of the broken symmetry for the MPEP, we have also 
followed the dynamics of the relaxational part of the es- 
cape paths, after they have crossed the y-axis separa- 
trix. The prehistory and relaxational probability distri- 
butions provide a complete history of the time evolution 
of large fluctuations to and from a given remote state. 
One can thus investigate experimentally detailed balance 
and time-symmetry (or the lack of them) |l^ . The posi- 
tions of the maxima of the measured relaxational distri- 
butions are compared with the corresponding theoretical 
trajectories in Fig. 2(b). A detailed analysis of the distri- 
butions will be given elsewhere. It can be seen from the 
figure that for a > ac the MPEP breaks time-reversal 
symmetry, i.e., the average growth and average decay of 
fluctuations |2^ traced out by the ridges of the corre- 
sponding distributions take place along trajectories that 
are asymmetric in time. That is, for a > ac the MPEP 
is not a time-reversed relaxational trajectory. 

The inset in Fig. 2(b) shows the distribution of points 
where the escape trajectories hit the y-axis separatrix 
(i.e., the exit location distribution). Its shape is nearly 
Gaussian, as expected from the saddle point approxima- 
tion of p^ . The maximum is situated near the saddle 
point clearly demonstrating that, in the limit of weak 
noise, exit occurs via the saddle point. 

The relationship between time-reversal symmetry- 
breaking for the MPEP when a > Uc, and symmetry- 
breaking generally for the system (0), is quite subtle. The 
system loses detailed balance and time reversal symmetry 
as soon as a > 1 and the drift field u becomes nongradi- 
ent. It is on account of a special symmetry of the system 
(reflection symmetry through the x-axis) that the MPEP 
can remain unchanged in this nongradient drift field up to 
the value ac — 4. Thus, for 1 < a < 4 the dynamics of 
the most probable fluctuational trajectories is a mirror- 
image of the relaxational dynamics only along the x-axis; 
everywhere else in the domain of attraction of (xg, 0) the 
outward optimal trajectories are not antiparallel to the 
inward relaxational trajectories, and the resulting closed 
loops enclose nonzero area pl] , |T^ . 

This prediction has been tested experimentally by trac- 
ing out optimal paths to/from specified remote states 
both on and off the x-axis, for 1 < a < ac- Some re- 
sults are shown in Fig. 3 for a — 3.5. It is evident that 
the ridges of the fluctuational (filled circles) and relax- 
ational (pluses) distributions follow closely the theoreti- 
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FIG. 3. Demonstration of local properties of detailed bal- 
ance and time-reversal symmetry of (|l|) for a = 3.5. (a) Mea- 
sured positions of the ridges of the fluctuational (filled circles) 
and relaxational (pluses) parts of the trajectories from (1,0) 
to (0.3,0.3), compared with theoretical predictions [ [Ll[|l2{ 
(curves), (b) The same for trajectories extending to the 
on-axis remote state (0.1,0). 



cal curves. For an ofF-axis remote state (Fig. 3(a)), they 
form closed loops of nonzero area, thus demonstrating 
the expected rotational flow of the probability current in 
a nonequilibrium system | pO[ . The corresponding ridges 
for an on-axis remote state (Fig. 3(b)) are antiparallel, 
indicating that symmetry is preserved along the x-axis. 

Our results verify the predicted bifurcation of the 
MPEP in (0) at a = ttc = 4, with a corresponding nonan- 
alyticity of the activation energy. We have demonstrated 
that, in the limit I? ^ 0, detailed balance and time- 
reversal symmetry can be considered as local properties 
along the MPEP of the system in the sense discussed 
above, and that the bifurcation phenomenon can be 
related to local time-reversal symmetry-breaking along 
the MPEP: results that may bear on two-dimensional 
stochastic ratchets pl| ] where symmetry plays an impor- 
tant role. Having thus demonstrated (see also [|l^) the 
reality of phenomena inferred from D — > optimal paths, 
we anticipate that other important theoretical 

predictions, e.g. "cycling" of the exit location distribu- 
tion ||2^, will also be physically realisable. 
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